`bsa` <-
function(wt,ht)signif(exp(-3.751) * ht^0.422 * wt^0.515,3)


